Effects of a quantum measurement on the electric conductivity: Application to graphene 
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We generalize the standard linear-response (Kubo) theory to obtain the conductivity of a system that is subject 
to a quantum measurement of the current. Our approach can be used to specifically elucidate how back-action 
inherent to quantum measurements affects electronic transport. To illustrate the utility of our general formalism, 
we calculate the frequency-dependent conductivity of graphene and discuss the effect of measurement-induced 
decoherence on its value in the dc limit. We are able to resolve an ambiguity related to the parametric dependence 
of the minimal conductivity. 
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The fact that measurements exert a back-action on the mea- 
sured object has attracted a lot of attention, 1-3 partly due to its 
relevance for the foundations of quantum physics, but also be- 
cause of implications for metrology 4 and the design of solid- 
state devices. 5 Fundamental considerations necessitate a dis- 
tinction between selective and nonselective descriptions of 
quantum measurements. 6 Selective descriptions use stochastic 
differential equations, 7 ' 8 or restricted path integrals, 9 and re- 
sult in conditional quantum dynamics when the measurement 
results are recorded. Some properties of selective measure- 
ments have been verified experimentally for areas as diverse 
as cavity QED 10 and superconducting phase qubits. 11 Nons- 
elective descriptions represent the evolution of the measured 
system irrespective of the measurement result. This descrip- 
tion takes into account all possible readouts, and the actual 
readout is assumed not to be known. 6,8 Quantum-mechanical 
back-action on the unsharply measured system causes loss of 
coherence between eigenstates of the measured quantity. In 
this work, we discuss measurement back-action theoretically 
within the nonselective framework. This approach makes it 
possible to determine how a macroscopic observable such as 
the conductivity of a system is affected when the current is de- 
tected in an unsharp-measurement scenario. We derive a gen- 
eralized Kubo formula where the measurement back-action 
provides a natural damping mechanism. We demonstrate 
the power of the developed general formalism by calculating 
the frequency-dependent conductivity <t(ui) of graphene^ a 
promising candidate for future micro- and nanoelectronics 13 
and also a low energy laboratory of relativistic physics 
Electronic -transport properties of graphene were analyzed in 
several previous studies using different methods; e.g., the 
Landauer-Biittiker formalism, 1617 the linear-response Kubo 
formula, 17 20 and the Boltzmann equation. 21 It was found that 
different parametric dependences of the dc conductivity can 
result from different limiting procedures applied to ordinary 
Kubo formulae. 20 Within our generalized Kubo formalism, 
we obtain physical conditions for when these results apply. 
Two regimes can be distinguished by comparison of the en- 
ergy scale HT that quantifies measurement-induced decoher- 
ence with the greater one among the thermal energy k^T and 
the chemical potential p (measured from the Dirac point). 



Weak back-action (HT < max{fcBT, A*}) results in Drude-type 
behavior <t(0) oc 1/T. In the opposite (strong-back-action) 
limit, a mixing of intra-band and inter-band contributions oc- 
curs that changes the parametric dependence of the dc con- 
ductivity such that <t(0) oc V. 

We employ the linear-response (Kubo) formalism 22 and di- 
vide the system's Hamiltonian into the part Hq, which gov- 
erns the free evolution, and SH, the perturbation associated 
with an external electric field E. For simplicity, we take the 
latter to be constant in space and assume the field to be applied 
between t — — oo and t = 0. The perturbation Hamiltonian 
is SH — — eE • f e ttJJt . In a nonselective description^! 3 - the 
dynamics of the density matrix, when the current is measured, 
is governed by a master equation with the back-action caused 
by a term of Lindblad form 23 



[H,p]- l\j,\j,p]] = £p-kSH,p], (1) 



-Ji[ H o,p] ~ b',/5]] with the current op- 
-[Ho,r] being our measured observable. 2 ! We 



dp 
~dt ~ 

where Cp 
erator j 

note that the measured observable is arbitrarily chosen to be 
the current, as it is relevant for the conductance calculation at 
hand. The main parameter of an unsharp quantum measure- 
ment device, as described by Eq. ([1), is the detection perfor- 
mance 7 = (Ai)" 1 (Aj)~ 2 , where At is the time resolution 
or, equivalently, the inverse bandwidth of the detector, and Aj 
the statistical error characterizing unsharp detection of the av- 
erage current. Within linear-response theory, we can linearize 
p = po + Sp, where po is the system's equilibrium density 
matrix. Keeping only linear terms in Eq. (fTJ, we get 



^ = CSp- ~[5H,p ] ~ ~L?\ d,Po]] > 



(2) 



assuming that the unsharp detection does not affect the equi- 
librium and using [Hq, po] = as well as [SH, Sp] ~ 0. In- 
troducing Ap = e~ ct Sp yields 



dAp 
dt 



-l[5H,po}-l\j,d,p }}) ■ (3) 



2 



Note that A/5 and 5p have the same value at t — 0, and both 
vanish at t = — oo. Integration yields 



= 0) = / dt e~ a [ - l -[5H,po] - -[?, [j,Po 



(4) 

The exponential factor in Eq. (0]i ensures convergence of 
the time integral, making it unnecessary to introduce the 
phenomenological adiabatic damping parameter employed in 
conventional linear-response theory, 22 Inserting (|4|i into the 
expectation value for the current density j and dividing by |E| 
yields the optical conductivity 



a fiU (cj) = [K^(t)e 1 ^ +K' liv {t)] dt , (5) 

J — OO 

with the kernels 

K^{t) = -ilr {fe-tQJ* ([ef v ,p ])} , (6) 



^W=Tr |e-^ (-1%, 



3^ Pol 



(7) 



Equations (O and (0 are the general result for the response 
to an unsharp measurement and could, in principle, be applied 
to any quantum system. As an instructive application, we 
now use Eq. (0 to calculate the frequency-dependent conduc- 
tivity of single-layer graphene within the continuum model 
for quasiparticles close to the K point in the Brillouin zone. 
In Eq. (Q]l, the Hamiltonian Hq could, in principle, include 
Coulomb-interactions, impurity scattering, density dependen- 
cies etc. In the present work, we just use the free-quasiparticle 



Hamiltonian for graphene in plane-wave representation, 

H (k) = hv(a x k x + (Jyky) , 



(8) 



where k x ,k y are the Cartesian components of wave vector 
k, <7, denote Pauli matrices acting in the sublattice-related 
pseudo-spin space, and v is the Fermi velocity, which has a 
value ~ 10 6 m/s. With the position operator f being the wave- 
vector gradient, the current operator is j M = ¥ [Hq (k) , ffj\ = 

f d ^j, ■ Single-particle eigenstates of clean graphene can 
be written as a direct product of a plane wave in configuration 
space with a spinor \n) = |k) <g) |cr)k- Here a labels the elec- 
tron and hole bands, respectively, and the spinor wave func- 
tion depends on wave vector k. The current operators with 
the Hamiltonian <[8J m the spinor space are j x = eva x , and 
j y = evayM From the definition of the equilibrium density 
matrix in the spinor space we find po|c)k = /(ftek,<r)|o')k, 
and Ho\c)-k = ^ek,(r|c r )k, where / is the Fermi-Dirac distri- 
bution function and ek.± = ±|k|. (For simplicity, the speed v 
has been absorbed into k.) Together with the anti-symmetry 
in momentum space, this implies K' = 0. The calculation 
of is straightforward, and using the Laplace transform to 
solve for the dynamics, we find 



d 2 k 



Res{^(k, 7 ,z)e zt }, (9) 



h J (2tt) 2 

where Res stands for the sum of residues of the integrand. It 
can be seen that K xy = K yx = 0. The remaining two conduc- 
tivities are identical, as a change of variables k x O k y in the 
expression for K xx yields K yy . The choice of conductivity 
measured will decide what kind of back-action will influence 
the system. Here we take j x , the current along the x-direction, 
which breaks the isotropy of the problem and we find 



K xx (k,T,z) 



k 2 x | k | (i6r 2 - 8zr + 4 | k | 2 +z 2 ) J2 a (- 



df(he) 



fc 2 (z-4r) 2 [/(fek,-)-/(ri6 k , + )] 



k | 3 [l6zr 2 - 8(2/c 2 + z 2 )r + z(4 I k | 2 +z 2 )] 



(10) 



r 



where the parameter T = r ye 2 v 2 /8 was introduced. The cubic 
factor in the denominator of Eq. ( fTOb has three roots Zj, which 
give the poles in Eq. ([9). 

For small V, the roots are to lowest order z\ = and 22,3 = 
±2i|k|. Using this and performing the time-integration in the 
limit r — > 0, the known intra- and inter-band contributions 
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8c osh(^)+co S h(^ r )' 



(12) 



to the conductivity of clean graphene are found. The scale fac- 
tor do = 4e 2 /h accounts for spin and pseudospin degeneracy. 

In the following, the conductivity is calculated numerically 
for finite values of T from Eq. (0 with Eqs. (O and ( fTOb . We 
use ksT as unit of energy. Figures Q] and [2] show the ac (opti- 
cal) conductivity, 27 whereas Fig.[3]shows the dc conductivity 
that is measured, e.g, in mesoscopic transport experiments. 

In Fig.Q] the conductivity is shown as function of frequency 
for different values of the coupling strength when the chem- 
ical potential remains fixed. For high frequencies, the con- 
ductivity saturates to the universal value 7r/8, indicated by 
a dashed line. The detailed shape of the crossover to satu- 
ration depends on T, with higher values pushing it to higher 
frequencies. In the limit of small T, measurement-induced de- 
coherence appears to simulate the effect of life-time broaden- 



3 



CO 










10 1 










10° 






fir , 

= 10 2 










\ ^<^<^' 




1CT 1 










io- 2 


io- 1 




io 2 




25 Ttu/HbT 



FIG. 1: Frequency dependence of the conductivity of clean single- 
layer graphene when the current is unsharply quantum measured. 
The chemical potential is fixed at fi/k B T = 1, and the value for 
measurement-induced decoherence assumed for each curve is indi- 
cated. The dashed line is at n/8. As the coupling to the measure- 
ment device is decreased, the curves more closely resemble the result 
found for clean graphene. 



ing due to inelastic scattering, 19,20 ' 26 but a closer look reveals 
that it is fundamentally different. Technically, the effect of 
r goes beyond merely broadening of distribution functions, it 
also moves the position of their peaks in energy, thus changing 
the resonance condition. For clean graphene, there is only a 
delta function peak for the mtra-band contribution, whereas 
here the intra- and inter-band contributions to the conductiv- 
ity become mixed. The existence of such a mixing has been 
inferred from recent experiments. 27 The dependence on the 
chemical potential is illustrated in Fig. [2] Generally, increas- 
ing the chemical potential shifts the frequency beyond which 
the conductivity attains its universal saturation to higher val- 
ues. This behavior is as expected theoretically 26 and observed 
in experiments. 27 At frequencies uj smaller than a few times 
the chemical potential, there is a significant departure from the 
universal conductance plateau. For fixed /i and k%T the sat- 
uration point in scattering models is independent of the scat- 
tering parameter, whereas in our work it strongly depends on 
the value of T. For /i > HT, the saturation occurs as in clean 
graphene, whereas the opposite case gives saturation for larger 
frequencies with increasing T, as seen in Fig.Q] 

When the effect of disorder is modeled conventionally by 
a life-time broadening due to inelastic scatteringr2i2£ a Drude 
peak is found for hio/k^T ~ 0, with a height inversely pro- 
portional to the inelastic-scattering rate. Since the master 
equation describing the effect of quantum measurements is 
formally identical to certain models of decoherence, we ex- 
pect that (t(0)/(7o ~ 1 /r. This turns out to be correct only for 
hT < max{/i, k^T}, as can be seen in Fig. 03. In the oppo- 
site limit, the parametric dependence of the dc conductivity is 
changed; it then grows linearly with T. A similar anomalous 
behavior was obtained in Ref. 20 by applying an unconven- 
tional limiting procedure within the Kubo formalism. Here 
we are able to readily identify the regimes where ordinary- 
Drude behavior or anomalous band-mixing behavior will be 
exhibited. Clearly, observing the latter should be easiest when 



FIG. 2: AC conductivity of single-layer graphene, for fixed 
hT/k B T = 1 and different values of chemical potential fi (mea- 
sured from the Dirac point). 



the chemical potential is at the Dirac point and the tempera- 
ture is low enough to satisfy KT > k^T, which corresponds 
to the minimal conductivity regime for graphene. Figure 3a 
shows how the nonmonotonic T-dependence of the conduc- 
tivity would be manifested in a typical transport experiment 
where the chemical potential (i.e., the density) is varied. 

The theory presented here is based on an unsharp measure- 
ment of the current density. To estimate the magnitude F of 
measurement-induced decoherence, we must consider the two 
situations most closely related to our result, optical conductiv- 
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FIG. 3: a) DC conductivity plotted as a function of chemical potential 
for different values of V. The conductivity is linear in the chemical 
potential for large values of fi/ksT. b) DC conductivity shown as 
a function of hV /k B T for different values of the chemical potential. 
For fi > /cbT and KT < \i, the conductivity is inversely proportional 
to the effective rate of decoherence, in analogy with the behavior 
expected from inelastic impurity scattering. In contrast, for KT > fi 
there is direct proportionality. This behavior was exhibited as long 
as /x > ksT, whereas curves for fj, < k%T practically coincide. 
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ity measurements, and mesoscopic transport measurements. 
For the latter case, it can be easily seen that e 2 u 2 /(Aj) 2 is the 
signal-to-noise ratio. Using typical values Av > 10 6 Hz for 
the bandwidth and AI / 1 « 10~ 3 for the signal-to-noise ra- 
tio, we find that values up to hT/k^T w 10 _1 — 10 2 can be 
achieved for T ~ 1 — 300 K. For measurements of the optical 
conductivity, the situation is more complex, as the measured 
signal is induced by the currents generated by the applied op- 
tical field. As a result, additional uncertainties such as geo- 
metrical factors and detector efficiencies become important, 
possibly bringing down the detection performance, but this 
can in principle be compensated by the increase in bandwidth 
offered by optical detectors, Av > 10 9 Hz. 

In conclusion, we derived a new Kubo formula to study the 
effect of measurement-induced back-action on the conductiv- 



ity. The back-action naturally introduces a source of damp- 
ing and thus makes the converged adiabaticity parameter fre- 
quently used in Kubo formula calculations superfluous. We 
applied this approach to calculate the electric conductivity of 
single-layer graphene. Mixing of the intra- and mter-band 
contributions to the dc conductivity strongly affect its para- 
metric dependence on the detector performance T. The regime 
of weak coupling to the measuring device models a standard 
Drude-type behavior, whereas in the opposite limit of strong 
back-action, we find that measuring a current in graphene will 
actually enhance the conductivity. 
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Royal Society of New Zealand. 
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